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ABSTRACT 



Context. The 'seeds' of the supermassive black holes (SMBHs) with masses of ~10 9 M Q observed already at z ~ 6 may have 
formed through the direct collapse of primordial gas in T viT > 10 4 K halos, whereby the gas must stay hot (~10 4 K) in order to avoid 
fragmentation. 

Aims. The interplay between magnetic fields, turbulence, and a UV radiation background during the gravitational collapse of primor- 
dial gas in a halo is explored, in particular the possibilities for avoiding fragmentation. 

Methods. Using an analytical one-zone model, the evolution of a cloud of primordial gas is followed from its initial cosmic expansion 
through turnaround, virialization and collapse up to a density of 10 7 cm 4 . 

Results. It was found that in halos with no significant turbulence, the critical UV background intensity (7,"') for keeping the gas hot 
is lowered by a factor ~10 for an initial comoving magnetic field So ~ 2nG as compared to the zero-field case, and lowered even 
more for stronger fields. In turbulent halos, /H" is found to be a factor ~10 lower compared to the zero-field-zero-turbulence case, 
and the stronger the turbulence (more massive halo and/or stronger turbulent heating) the lower y™'. 

Conclusions. The reduction in 7™' is particularly important, since it exponentially increases the number of halos exposed to a 
supercritical radiation background. 
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1. Introduction 

Several very bright quasars have been detected already at z > 6. 
This suggests that some SMBHs with masses of =* 10 9 M q al- 
ready existed when the Universe was less than 1 Gyr old (tFanl 
2006). Explaining how such massive SMBHs could have as- 
sembled so soon after the Big Bang presents quite a challenge. 
The main questions concern how and when the 'seeds' of these 
SMBHs formed and how their subsequent growth proceeded. 

Several pathways leading to the formation of seed black 
holes have already been proposed. One group of scenarios sug- 
gests that seed black holes formed via the direct collapse of 
metal-free (or very metal -poor) gas in halos with T vlr > 10 4 K, at 
redshifts ~5 - 1 0, resulting in seed black holes with M ~ 10 4 - 
10 5 Mr, (see e.g iBromm & Loebl200llKoushiappas et al/2 004; 
Begelman et al]l2006l: [Lodato & Nataraianl l2006: Sp aans & Silkl 



2006t [Schleicher et al.ll2010bl) . For efficient gas collapse to oc- 
cur, fragmentation must be suppressed. This is possible if the 
gas in the halo is kept hot (and thus the Jeans mass high). Hence, 
the formation of H 2 must be inhibited so cooling can occur only 
through atomic hydrogen, as otherwise H 2 cooling will lower 
temperatures to ~200 K. 

Several mechanisms have been suggested to suppress H 2 
cooling and thus prevent fragmentation. The main one of these 
mechanisms requires a critical level of Lyman-Werner radiation 
to photo-dissociate H 2 and keep its abundance very low. The 
critical intensity needed to suppress H 2 in the massive halos 
where direct gas collapse can occur is large compared to the ex- 
pected cosmic UV background at the relevant redshifts. How- 



ever, the distribution has a long bright-end tail and halos irradi- 
ated b y supercritical inten sities would be a small subset of all 
halos (iDiikstra et al.ll2008l) . 

Another mechanism proposes that the dissipation of a suf- 
ficiently strong magnetic field can heat the gas in the halo to 
~10 4 K, which causes H 7 to be destroyed by collisional dissoci- 
ation (ISefhi et al.ll2010h . 

A variety of mechanisms exist for generating magnetic fields 
early in the Uni verse, both before an d after recombination (for a 
review, see e.g. IWidrow et al.ll2012l) . and also for amplifying an 
existing field. In the case of a collapsing halo, the most important 
ones are gravitational compression, the small-scale turbulent dy- 
namo, the large-scale dynamo in protostellar and galactic disks, 
and the magneto-rotational instability (MRI). 

Gravitational compression increases the magnetic field as 
B ex p® when the field is coupled to the gas. For spherically sym- 
metric collapse, a = 2/3, but if the collapse proceeds preferen- 
tially along one axis, the scaling is closer to a - 1 /2. In realistic 
cases , often intermediate values are found (e.g. ISchleicher et all 
l2rjQ^lHoaAetd]l20Tl . 

Non-helical turbulent flows can act as small-scal e dynamos, 
which produce disordered, random magnetic fields (Kazantsev 
1 19681) . The magnetic field amplification results from the ran- 
dom stretching and folding of the field lines by the turbulent ran- 
dom flow. During gravitational collapse, turbulence is generated 
by the release of gravitational energy and the infall of accreted 
gas on the inner, self-gravitating core. This means that, in the 
context of star and galaxy formation, a strong tangled magnetic 
field may be generated already during the collapse phase by the 
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small-scale dynamo (ISchleicher et al.ll2010ah . For the formation 
of seed black holes, it implies that the existence of an accretion 
disk may cause the magnetic field to be further amplified (by a 
large-scale dynamo and/or the MRI) which provides additional 
stability and hence reduces fragmentation. 

2. The model 

The evolution of a cloud of primordial gas is followed from ini- 
tial expansion to a high-density core, using a one-zone model, 
in which the physical variables involved are regarded as those at 
the cloud center. The model uses standard cosmology, with c os- 
mological parameters given by WMAP7 dLarson et al.ll201 lb . 

2. 1 . Density evolution 

The spherical collapse model for a top-hat overdensity is used 
to compute the matter density. At the moment of turnaround, 
the gas decouples from the dark matter and becomes self- 
gravitating. Any effects due to rotation are neglected for sim- 
plicity. The baryonic matter collapse is expected to proceed 
like t he Larson-Penston similarity solution (for the isothermal 
case; iLarsonl 1969t Penston| [ l9 69), as generalized to polytropic 
cases bv lYahilld 1983b . According to this solution, the cloud con- 
sists of two parts: a central core region and an envelope. The 
central core region has a flat density distribution, whereas the 
density in the envelope decreases outwards. The size of the 
central flat region is roughly given by the local Jeans length, 
/lj = c s yjn/ (Gp m ), with c s - sJykftT / (/jniu) the sound speed 
in the central region and p m the total matter density in the cen- 
tral region. The collapse in the core proceeds approximately 
at the free-fall rate, although additional heat input, for exam- 
ple due to magnetic energy dissipation, may delay gravitational 
collapse. The mean baryonic density evolution in the central 
part is described by p\ - pt,/fff where % is the free-fall time: 
t«= V3tt/ (32Gp m ). 



2.2. Chemical network 

The species that are included in the chemical network of this 
model are H, H + , H~, H 2 , H 2 , and e~. HD or other molecules 
involving deuterium are not included. Reactions with He are 
not taken into account, but He is considered in the calculation 
of the mean molecular mass. The He mass fraction is taken to 
be ~0.248 (corresponding to an abundance xu e ~ 0.0825) and 
stays constant throughout the time integration. The fractional 
abundances of H, H, and e~ are explicitly followed during the 
integration. 

The evolution of the fractional abundan ce of electrons , x e , 
is given by the following equation (see e.g. lPeebleslll993l for 
further details): 



where 



dx e 
~df 



B e x m exp 



hv„ 



a e x e n H 



C + y e (T)x m x e n H , (1) 



ks T cm b 

The first term represents the recombination and photo-ionization 
of the primordial plasma, the second term is the collisional re- 
combination term and the third term represents collisional ion- 
ization (H + e~ — > H + + 2 e"). 

The evolution of the fractional abundance of m olecular hy- 
drogen, x^, is given by the following equation (IShang et al.1 
l20ialSefliietalJl2010h : 

dx H , 



kgkioxmnn 



kmxmnu +k y + {k 13 + fc 2 i) *e«H 

+ k l9 x e n H + k 2 Qx m nn + k 25 

^des = k l5 X m + k n X p + k l& X e + £ 2 8/sh/«H- 



(3) 

(4) 



Here, k m is the net rate of formation of H 2 through the H chan- 
nel, £des is the net destruction rate of H 2 , and k y is the destruction 
rate of H by CMB photons. The reaction rates can be found in 
appendix A of lShang et al.l (1201 Ol) . numbered as above. 
For column densities that are large enough, H 2 can shield itself 
from radiation in the Lyman- Werner bands. The self-shielding 
facto r / s h is given by Draine & Bertoldi (1996). The updated / s j, 
from I Wolcott-Green et al.l (1201 ll) results in decreased shielding 
and may thus enhance our results. 



2.2.1. Radiation background 

A sufficiently intense UV radiation background can either di- 
rectly photo-dissociate H 2 , or photo-dissociate the intermediary 
H~. The relevant criterion for suppressing H 2 formation is that 
the photo-dissociation timescale is shorter than the formation 
timescale, which results in 7 cnt oc p. The intensity is written 



At 



= k m x e x Hl n H - kfe s x H ,n H , 



(2) 



as 721, which denotes the specific intensity just below 13.6 eV, 
in the units of 10 ergcm~ 2 sr ~' s Hz -1 . Here, two differ- 
ent UV spectra are considered. They are both Planck spectra 
with a blackbody temperature of either T* = 10 4 or 10 5 K (T4 
and T5, respectively). The T4 spectrum is meant to approxi- 
mate the mean spectrum of a normal stellar population (Pop II), 
whereas the T5 spectrum is closer to the harder spectrum ex- 

S iected to be emitted by t he first generation of stars (Pop II I) 
Tumlinson & Shullll2000t iBromm et al1l200lb ISchaereril2002l) . 

2.3. Magnetic field evolution 

Gravitational compression and the small-scale dynamo can am- 
plify the magnetic field strength B, while ambipolar diffusion 
will decrease B. If the flux-freezing condition applies, the mag- 
netic field depends on the density as B oc p?, where a lies in 
the range 2/3 - 1/2. Hence, the magnetic field strength will 
increase during gravitational collapse. It is assumed that gas 
is continually falling in, so the turbulence generated by accre- 
tion will not decay but is constantly replenished. However, 
this depends on the ambient gas reservoir and may in reality be 
more complicated. It has been shown that the injection scale of 
such accretion-driven turbulence i s close to the size of the sys- 
tem under considerati on (thus, As) (IKlessen & HennebellefeOlQ : 
iFederrath et al.ll201lb . The turbulent velocity on the injection 
scale is expected to be comparable to the velocity of the in- 
falling gas, and for a roughly isothermal density profile, the free- 
fall velocity is independent of radius. Hence, while the injec- 
tion scale changes during the collapse, the injected velocity v\ n 
sta ys the same and is approximately equal to the virial v eloc- 
ity dGreif et al.ll2008t IWise & Abeill2007l IwTse et alJl2008h . On 
scales smaller than the injection scale, the turbulent velocity is 
expected to scale as v oc f, with/? =1/3 for Kolmogorov turbu- 
lence (incompressible gas) and B — 1/2 for Burgers turbulence 
(strongly compressed gas). 

The magnetic field on a scale I typically grows exponentially 
on the eddy turnover time f ec j = l/v, where v is the turbulent 
velocity on the scale I. However, the magnetic field will satu- 
rate when the magnetic energy corresponds to a fraction Rm~} of 
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the kinetic energy. Rm a is the critical magnetic Reynolds num- 
ber; if Rm > Rm cr , the magnetic field grows through dynamo 
action. The maximum magnetic field strength is thus given by 
B max = yj4npbV l /Rm CT dSubramanian & Barrowlll998b . How- 
ever, the exact value of the saturation field strength is still some- 
what uncertain. Once B > B max on a certain scale, it is no longer 
amplified by the small-scale dynamo; however it is still ampli- 
fied by gravitational compression (ex p a ). It can thus in principle 
increase above the saturation level (which only increases ex p 1/2 ), 
but in this case it is subject to turbulent decay. On a given scale, 
this decay will probably also occur on f^j. As a result, the mag- 
netic field strength tends to stay close to B max on that scale. 

The most important contribution to the total magnetic energy 
comes from the integral scale, the scale on which the magnetic 
field is largest. ISchleicher et al.l (1201 Pal) have shown that in an 
atomic cooling halo, the integral scale increases very rapidly to 
the maximum scale on which the magnetic field is coherent after 
the start of the simulation. For this reason, here only the evo- 
lution of the magnetic field at this scale of maximal coherence 
is followed. Since the magnetic field is distorted by the gravi- 
tational collapse, the largest possible coherence length is always 
smaller than the Jeans length by some factor f&\ its precise value 
is uncertain. We adopted a fiducial value of 0.1; changing this 
by a factor of a few does not significantly affect the results. 

Ambipolar diffusion (AD) is important in a mostly neutral 
medium where a tangled magnetic field is present. The AD heat- 
ing rate can be estimated as Lad ~ ?1adB 2 / (47r/g) dShang et al.l 
120021; Sch leicher et al. 2008 ), where t jad is the AD resistivity 
(iPinto et al.ll2008t iPinto & Gallill2008h and / B is the coherence 
length of the magnetic field, which is calculated as the minimum 
of the Alven damping scale and the integral scale. 

The evolution of the magnetic field energy E B = B l j%n is 
then calculated as: 



T4; J„. = 10 



T5; J„, = 10* 



d£ B 

At 



2aE B lnp' b - L AD 
2£ B (lnfl) d 



/dynamo 



Z > Zvii, 
Zvir > Z. 



2.4. Thermal evolution 

The evolution of the temperature, T, is given b 



equation (see e.g. lPeeble 



iperature, 
il99l lS 



efhi et al.ll20O 



(5) 



the following 



dT 
~d7 



y-i 
Pb 



r^ + _ h -L c ) 

at Kb 



+ k c x e (T cmb - T) . (6) 



The first term in this equation is the adiabatic heating/cooling 
rate due to collapse/expansion, the second term incorporates 
various other heating (magnetic energy dissipation by AD and 
turbulent energy dissipation, as discussed below) and cooling 
(through H 2 and H) volume rates, and the third term represents 
Compton heating/cooling. 

Part of the turbulent energy will go to driving the small- 
scale dynamo, and part of it will be transferred from large 
scales to smaller ones in a cascade process, until it is dissi- 
pated by viscosity at small enough scales. The rate (per unit 
mass) at which energy i s injected into the system is e m - 
Einlmlt ei (A } ) = vr ;/(2/tT) (IShulll992l) . The volume heating rate 
from the accretion-driven turbulence (ADT) can then be esti- 
mated as Ladt = ftPb^in, where / t is the fraction of the injected 
energy that is dissipated. Its fiducial value is taken to be 0.1, 
changing this by a factor of a few does not significantly affect 
the results. 




10 10 10 
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Fig. 1. Gas temperature as a function of density for a T4 background 
with J21 = 10 (left) and a T5 background with J21 = 10 4 (right), for 
several different initial magnetic field strengths (no turbulence). The 
dotted vertical line indicates virialization. 



The model is initialized at a redshift of 800 (zinit)- The dissipa- 
tion of magnetic fields into the IGM after recombination can sig- 
nificantly influence its temperature and ionization, so this early 
start provides the proper initial conditions. We follow the evo- 
lution of a 10 9 M Q halo that virializes at z = 10. The radiation 
background is switched on at turnaround, z ~ 15; an earlier or 
slightly later turn-on is also possible, but does not change the 
results significantly. The integration is stopped when the density 
reaches n\, ~ 10 7 cirT 3 . The model is not equipped to deal with 
higher densities; this would require the inclusion of additional 
physical processes, e.g. three-body interactions which increase 
the H 2 formation rate. 



3. Results 

After the radiation background is turned on, H 2 is destroyed 
rapidly. It cannot self-shield as the density is too low; this hap- 
pens even for a low intensity of /21 = 1- The cooling is domi- 
nated by H, so the temperature is high and the ionized fraction 
becomes elevated, which in turn aids the formation of H 2 . For 
intensities that are low enough, H 2 succeeds in reforming, and 
becomes the dominant coolant. The density at which this occurs 
depends on the radiation intensity, and is higher for higher in- 
tensities. However, at a certain intensity H 2 cannot reform fast 
enough to become an important coolant. If turbulence is unim- 
portant and the initial magnetic field is <0.01 nG (virtually iden- 
tical to the zero-field case), this critical intensity is found to be 
10 < J 2 \ < 10 2 for a T4 background and 10 4 < J 2 \ < 10 5 for a 
T5 background. Similar results are found for Bo = 1 nG, always 
assuming a = 2/3. However, if the magnetic field is increased 
to 2 nG, H 2 never becomes an important coolant for a T4 back- 
ground with 721 = 10, and for /21 = 1 an instability occurs at 
high densities. Here, the AD heating becomes too strong to be 
compensated by the H 2 cooling and the temperature suddenly 
increases, as much of the H 2 is destroyed by collisional disso- 
ciation and H cooling is still very strongly suppressed below 
~8 x 10 3 K. The gas stays hot afterwards, with cooling domi- 
nated by H. For a T5 background with the same By, a similar 
instability occurs for J21 = 10 4 , while with Bo - 3nG the H 2 
fraction is never large enough for significant cooling at this in- 
tensity. Thus, when turbulence is not important, a larger initial 
magnetic field decreases the critical intensity required to keep 
the gas in the halo hot, as illustrated in Fig.Q] 

During collapse, turbulence quickly brings the magnetic field 
strength to B max through the small-scale dynamo, either by am- 
plifying smaller fields or by draining energy from larger fields. 

1 /2 

fi max only increases as <x n . , which means that heating from 
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Fig. 2. Gas temperature, electron fraction, and H 2 fraction as a func- 
tion of density for a turbulent halo with Bq = 1 nG, a T4 (left) and a T5 
(right) background, and different 7 2 i as indicated in the top panels. The 
dotted vertical line indicates virialization. 



turbulent dissipation always dominates over AD heating. H 2 
cooling grows steeper with density than turbulent heating, and 
thus the gas always cools through H 2 when there is no radiation 
background present. Since heating is dominated by turbulence, 
halos with different Bo converge to approx. the same evolution- 
ary track; the turbulence has a moderating effect. 

The effects of different radiation intensities of a T4 and T5 
background on the temperature, electron fraction and H 2 frac- 
tion are shown in Fig. [2] for a turbulent halo with Bo = 1 nG. 
However, the critical intensity does not depend on Bo when tur- 
bulent effects are important. The critical intensity is found to be 
1 < J 2 i < 10 for a T4 background and 10 3 < J 2 \ < 10 4 for a T5 
background. 

4. Discussion 

The effects of magnetic fields and turbulence on the critical in- 
tensity of a UV radiation background were examined. Note 
that these results only hold in a zero/very low metallicity en- 
vironment, as metals and dust are much more efficient coolants 
than H 2 . For a halo not significantly influenced by turbulence 
or magnetic fields, the critical intensity was found to be 10 1 < 
J c 2 f < 10 2 for a T4 background, and 10 4 < J%* < 10 5 for a 
T5 background. Th ese limits are consistent with those found by 
IShang et all (1201 Oh and lowe r by a factor ~ 10 than previo us esti- 
mates by e.g. lOmukail d200lh and lBromm & Loebl (J2003I) . likely 
due to the different H 2 dissociation rates used. For Bo = 1 nG, 
these limits do not change; however, when the field is increased 
to ~2 nG, the critical intensity required to suppress H, cooling 
is lowered by a factor ~10, and the stronger the field, the lower 
7™'. Such magnetic fields alone do not give rise to sufficient AD 
heating to overcome H 2 cooling, but in combination with a radi- 
ation background they do influence the ability of the gas to stay 
hot. Note that the amount of AD heating depends on the scaling 
of B with p; it is therefore important to obtain a correct model 
for this relationship. The current u pper limit on the primordial 
magnetic field is ~1 nG comoving (ISchleicher & Miniatill201 it 



iTrivedi et al.ll2.Q12h . so a 2 nG field would be reached by the ~2<x 
upward fluctuations. 

In a 10 9 M turbulent halo, the critical intensity was found to 
be 1 < J c 2 f < 10 1 for a T4 background, and 10 3 < /™ l < 10 4 for 
a T5 background; these are a factor ~10 lower than for a halo not 
influenced by turbulence or magnetic fields. Since this is due to 
the turbulent heating in such halos, larger halos and/or halos with 
stronger turbulent heating will have an even lower /£?*. Note 
that the results for the non-turbulent halos are independent of 
halo mass. Interestingly, in turbulent halos with M > 10 11 M 
(also depending on the strength of the turbulence), the turbulent 
heating alone is strong enough to keep the gas hot, without any 
UV background. 

The fact that the values of 7™ 1 that have been found here are 
smaller than previous estimates is quite important. The mean 
cosmic UV background is expected to be around J 2 f ~ 40, 
which is smaller than the critical intensities, especially in the 
case of a T5 background. However, the background will fluc- 
tuate spatially, and thus a fraction of all halos will be irradiated 
by a supercritical intensity. Then, the lower the critical intensity, 
the larger the fraction of halos which are suitable candidates for 
direct S MBH formation; e.g. according to the distribution pro- 
posed bv lDiikstra et al.l(l2008l) . a decrease in 7™' from 10 4 to 10 3 
means an increase in the fraction of irradiated T w \ r « 10 4 K halos 
from negligibly small (< 10~ 8 ) to ~10~ 6 . For /™'~10 2 , the halo 
fraction even increases to ~10~ 3 . With a sufficiently low 7™ 1 , 
one could argue that this mechanism provided many, if not all, 
seeds for the SMBHs observed in galaxies today. 
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